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Transition metal-phtahalocyanine(Pc) compound, TPP[Fe(Pc)(CN)2]2, which is one of 
molecular conductors of charge transfer type with 3/4- filled conduction band consisting of 
LUMO of Pc molecules, shows characteristic features in transport and magnetic properties re- 
sulting from localized magnetic moments 5 = 1/2 associated with Fe +3 atoms. We construct an 
effective tight-binding model of this system and study the mechanism of exchange interaction, 
J, between d and tt electrons based on both second order perturbation of transfer integrals 
between d and n orbitals and numerical diagonalization. It is found that there is no hybridiza- 
tion between d-orbitals and LUMO of 7r-orbitals and then super-exchange interaction in the 
Anderson model does not exist. Instead, processes associated with Hund's rule both on d and 
7r orbitals, which may be called "the double Hund's exchange mechanism", turn out to play 
important roles and the sign of resultant J can be either ferromagnetic or antiferromagnetic 
depending on model parameters because of competition among various processes. By taking 
account of magnetic anisotropy due to spin-orbit interactions and comparing with experimental 
results, it is indicated that J is antiferromagnetic and of the order of 100K. 
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1. Introduction 

Phthalocyanine (Pc) compounds with transition metal 
(TM) arc organic materials with a variety of geometry 
together with the interplay between conduction electrons 
and magnetic moments of TMs. They show many inter- 
esting phenomena such as charge ordering, giant magne- 
toresistance, Kondo effect, ferrimagnetism, and so on. 

Recently, it has been reported that one of the TM-Pc 
systems, TPP[Fe(Pc)(CN) 2 ] 2 , shows large negative mag- 
netoresistance. 1-3 Because the valence values are TPP 1+ , 
CN 1 -, Pc 3 / 2 " and Fe 3+ (3d 5 and the low spin state), this 
compound has a one-dimensional conduction band of 3/4 
filling constructed from a lowest unoccupied molecular 
orbital (LUMO) of Pc, and has local moments of S = 1/2 
coming from the d orbitals of Fe. On the other hand, 
a similar compound, TPP[Co(Pc)(CN) 2 ]2, shows large 
positive magnetoresistance. 4 This system has a similar 
one-dimensional conduction band as in Fe-Pc compound, 
while it does not have the local moment on Co (S — 0). It 
is also known that TPP[Fe .o7Coo.93(Pc)(CN) 2 ]2 shows 
negative magnetoresistance as in Fe-Pc. 5 From these re- 
sults, it is natural to think that the origin of the negative 
magnetoresistance is due to the interaction between the 
conduction electrons (tt orbitals on Pc) and the local mo- 
ment on Fe. We call this 7r-d interaction. 

Assuming that the ir-d interaction is antiferromagnetic 
(AFM), the origin of negative magnetoresistance was 
discussed on the basis of an extended Kondo model. 6 
On the other hand, recently it was claimed that the 
7r-d interaction is ferromagnetic (FM) on the basis of 
quantum chemical calculation. 7 On the basis of the fer- 
romagnetic 7r-d interaction, the electronic state of the 



extended Kondo model has been discussed. 8,9 Then, it 
is suggested that the origin of the negative magnetore- 
sistance is double exchange interaction between Fe and 
Fe. 10 However, the negative magnetoresistance has been 
observed in TPP[Feo.o7Coo.93(Pc)(CN) 2 ]2 as introduced 
above, and the result suggests that the origin of the nega- 
tive magnetoresistance is not the double exchange mech- 
anism, because the mean distance between Fe and Fe is 
about 70 A. 

Based on these situations, it is important to study the 
7r-d interaction microscopically. Although the amplitude 
and its sign of the ir-d interaction was evaluated by the 
quantum chemical calculations, it is suspected that the 
quantum chemical calculations only take account of the 
low energy states. Thus, we think that the sign and the 
amplitude of the interaction is not decisive in real com- 
pounds. The mechanism of 7r-d interaction itself also has 
been unclear in the quantum chemical calculation. 

In this paper, we study the electronic states of 
[Fe(Pc)(CN) 2 ] in detail, and clarify the mechanism of 
7r-d interaction. First, by analyzing the electronic state 
of Pc, we find that there is no hybridization between d- 
orbitals and LUMO of 7r-electrons. This means that there 
is no super-exchange interaction between d electrons and 
7T electrons. Instead, the mechanisms of exchange inter- 
action are "generalized" Kanamori-Goodenough mecha- 
nisms. In particular, we find that a second order pro- 
cess which uses two kinds of Hund's rule couplings both 
on the d and tt electrons plays an important role. This 
mechanism gives an antiferromagnetic exchange interac- 
tion even if it is not a super-exchange interaction. We 
call this "double Hund's exchange" mechanism. 
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We find that the sign of the ir-d exchange interac- 
tion can be both ferromagnetic and antiferromagnetic 
depending on the model parameters because of the pres- 
ence of the competing processes. However, we think that 
the sign of the n-d interaction can be determined by a 
comparison with experiments. When we discuss the com- 
parison with experiments, it is necessary to take account 
of the anisotropy of the 7r-d interaction. For this purpose, 
we consider a material, PNP[FePc(CN) 2 ], which has the 
same unit FePc(CN) 2 as in TPP[FePc(CN) 2 ] 2 . 2 The unit 
of FePc(CN) 2 has a local magnetic moment of S = 1/2 of 
Fe atom, while it has no tt conduction electrons, because 
the valence is PNP 1+ , CN 1 ", Pc 2 ~ and Fe 3+ . In this 
case, we can study the electronic state purely of d elec- 
trons on Fe site without tt electrons. Therefore we discuss 
the electronic state of d electrons in this PNP compound 
by analyzing the anisotropy of the g values observed in 
EPR experiment. On the basis of these results, we derive 
the effective model of 7r-d interaction, and discuss its sign 
by comparing with susceptibility obtained in the dilute 
system of Fe atoms, TPP[Coo.9 3 Feo.o7Pc(CN) 2 ] 2 . 

This paper is organized as follows. In §2, the electronic 
state of Pc is discussed in the tight binding approxima- 
tion based on tt orbitals on atoms in Pc. We show that 
the hybridization between LUMO and d orbitals will var- 
nish in this system. In §3, on the basis of the electronic 
states of tt orbitals clarified in §2, the exchange interac- 
tion between d electrons and tt electrons is derived and 
the amplitude of this interaction is estimated. In §4, to 
consider the anisotropy of 7r-d interaction, the electronic 
state of d electrons is studied by analyzing the g factor 
in PNP[FePc(CN) 2 ] on the basis of a simple model in- 
cluding a crystalline electric field (CEF) and a spin-orbit 
coupling. In §5, the effective model of 7r-d interaction is 
derived by using the results of §3 and §4, and the sign of 
the TT-d interaction is discussed by the comparison with 
experimental results. 

2. Electronic state of Phthalocyanine 

Figure 1 (a) shows the structure of metal- free phthalo- 
cyanine (Pc). Here Nj (i=l ~ 8) represents the i-th nitro- 
gen (N) site, and Ai ~ Ds are carbon (C) sites, respec- 
tively. The electronic state of Pc has been understood in 
detail by both experiments (XAS and XPS) and quan- 
tum chemical calculations. 11-16 

In the following, we show that these experimental 
and calculational results can be semiquantitatively re- 
produced by the tight binding approximation using tt 
orbitals. As defined in Fig. 1(a), the transfer integrals 
i(black line), and <i(blue line) are the transfer integrals 
between C atoms in hexagons, those between C atoms in 
pentagons, respectively. i 2 (red line), and £3 (green line) 
are those between C and N. We diagonalize this tight- 
binding model and adjust the parameters, t ~ £3 so as 
to reproduce the results of the quantum chemical calcu- 
lations. 

Figure 1(b) shows the energy levels of the tt orbitals 
when we choose t\ — 0.8i, t% = l.lt, t% = 1.2t, and 
An = —0.9t. This result is consistent with the quantum 
chemical calculations quite well. Here, An is the one- 
body level of N measured from that of C. 



As shown in Fig. 1(b), the gap between LUMO and 
LUMO+1, +2 is of the order of OM. Since t ~ 3eV, this 
gap is estimated as 1.5 eV, which is in good agreement 
with the quantum chemical calculations (1.8eV for i^ 2 Pc 
molecule, and 1.6 eV for LiPc). 12,16 It is noted that the 
energy difference between the highest occupied molecular 
orbital (HOMO) and HOMO+7 is tiny ( O.lt ~ 0.3eV). 
Indeed this difference is estimated as 0.14 eV for LiPc in 
the quantum chemical calculations. 16 
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Fig. 1. (Color online) (a) Structure of metal-free phthalocyanine. 
Ni (i=l ~ 8) is i-th nitrogen (N) site, and Ai ~ Dg are carbon 
(C) sites. t(black line), and <i(blue line) are the transfer inte- 
grals between C atoms in hexagons, those between C atoms in 
pentagons, respectively. t2(red line), and 43(green line) are those 
between C and N. (b) Energy level scheme of tt orbitals. We also 
introduce Ajv which is the one-body level of N measured from 
that of C. 



The wave function of the i-th tt orbital is expressed as 
follows: 



i'=l~8 



(1) 



where 1 



and i>.\ are the wave functions of the tt orbitals 



of j-th N and C, respectively. The coefficient, c*, is the 
amplitude of the wave function of j-th N or C atom in the 
i-th tt orbital. Some of them (c l N ., j =1,3,5,7) are shown 
in Table I. These coefficients of N sites on LUMO are 
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Table I. The amplitude of coefficients of re orbitals for the four 
central N atoms, i.e, Nj (j=l,3,5,7) sites. 



zero, while those on LUMO+1, LUMO+2, and HOMO, 
HOMO+1 and HOMO+2 are finite. As discussed later, 
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these coefficients are important to discuss the hybridiza- 
tion between the d orbitals in TM and 7r orbitais, and 
we find that there are no hybridization between d or- 
bital and LUMO. In addition, we neglect HOMO+3 ~ 
HOMO+7 in Table I, because the hybridizations between 
these orbitals and d orbitals are much smaller than the 
hybridization between LUMO+1 and d orbital as shown 
in §3. 

The schematic pictures of wave functions are illus- 
trated in Fig. 2. The size of circle on each atom cor- 
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Fig. 2. (Color online) Schematic pictures of wave functions 
of LUMO+l(<t>i), LUMO+2(<t> 2 ), LUMO(* 3 ), HOMO($ 4 ), 
HOMO+l($ 5 ), and HOMO+2($ 6 ), respectively. 



responds to the amplitude of the coefficient of the wave 
function, and the blue (black) circle indicates positive 
(negative) sign. The composition of these wave functions 
is consistent with the quantum chemical calculations. 

3. Mechanism and estimation of the exchange 
interaction between the d electron in Fe and 
tt electron in Pc 

In this section, we discuss the 7r-d interaction in 
TPP[FcPc(CN)2]2 using the wave functions of 7r-orbitals 
in Pc obtained in the previous section. Although in 
TPP[FcPc(CN) 2 ] 2 the nominal values are Fe 3+ and 
Pc~ 3 / 2 , we discuss the simple case of Pc~ and Fc 3+ , i.c, 
there is one electron in LUMO of Fig. 1 and five electrons 
in 3d orbitals of Fe atom (see Fig. 3). As a microscopic 
model, we consider the following cluster model Hamilto- 
nian constructed from the d orbitals of Fe and tt orbitals 
of Pc; 



H = H t + H H 



H 



Pc- 



(2) 



where H t contains the transfer integrals between the d 
and 7r orbitals, and Hd and Hp c arc Hamiltonians for the 
d and tt orbitals, respectively. 
First, H t in eq. (2) is given as 



6 

e 

»'=! J = l(j¥4) 



(3) 



where tij is the transfer integral between d^ orbital and 
&j orbital, where we define d\ a = d yzcri d 2(T = d zxa , 
and d 3a = d xya , and $i = LUMO+1, $ 2 = LUMO+2, 
$ 3 = LUMO, $ 5 = HOMO+1, and $ 6 = HOMO+2 
in Fig. 2, respectively. We have neglected the other tt 
orbitals, because their energy levels are much higher (or 
lower) than that of LUMO, and the hybridizations with 
d orbitals are much smaller than those of LUMO+1. 

There is no hybridization between LUMO (conduc- 
tion band) and d orbitals (localized spins). Therefore, 
the mechanism of ir-d interaction is different from the 
mechanism based on the super-exchange interaction as 
in the impurity Anderson model which uses hybridiza- 
tions. 17 Thus, we consider different mechanisms for the 
7r-d interaction in the Fe-Pc complex. 

In eq. (3), we estimate the transfer integrals as —tn = 
t 2 2 = Q.38t dN , ti2 = <2i = 0.2At dN , and ti 5 = t 26 = 
—0.4At d N, and — ti6 = £25 = 0.36tdN- Here, tdN is the 
transfer integral between d zx (or d yz ) orbital and p z (tt) 
orbital of the nearest-neighbor N atoms, and the coeffi- 
cients of tij have been determined from the amplitude 
of the wave functions of 7r-orbital c* in eq. (1) and from 
the symmetry of the t 2ff orbitals. It is important in the 
following considerations that d xya (da a ) orbital does not 
hybridize with any 7r orbitals due to the symmetry re- 
quirement: i.e., t^j — 0. Non zero transfer integrals are 
shown in Fig. 3 by the green dashed lines. 

Next, Hd in eq. (2) is given as follows: 



i— 1~3,<7 
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(4) 



(5) 



(6) 



Jd 
' 2 



E 

ij = l-3(i/j) 



i^djtd^ + h.c), (7) 



where Af 's are the one-body levels of the d orbitals as in 
Fig. 3. Ud, U' d and Jd represent the intra- and inter orbital 
Coulomb interaction, and the FM exchange interaction 
corresponding to the Hund's rule coupling. Here, we set 
these parameters as Ud — 5.0tdN, and Af = A d (i = 1 ~ 
3): i.e., the CEF and the spin-orbit coupling are neglected 
because the energy scale of the CEF and the spin-orbit 
coupling is much smaller than the Coulomb interaction 
between d orbitals. We will consider the effect of these 
terms to discuss the anisotropy of the tt- d interaction in 
§4 and §5. 

Finally, Hp c in eq. (2) is given as 



H 



Pc — 



J 



M 



A f"i-- E ~Y n J- n ^ (8) 



j=l,2,5,6,CT 



J 



.7=1,2,3,6 



M 
J_ 

2 



E 4-(^t*3t*^+h.c)»(9) 



where A*^ are the one-body levels of $1, $ 2 , and $3 
measured from $5 or $6, respectively, for Pc. The level 
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scheme of Aj 's is shown in Fig. 3 where we take Af 

i A/ — 



A/ 



Af = A L+ i = 3.5cV and Af = A L = 2.0eV because 
A L+ i ~ 1.18* - 3.54eV and A L ~ 0.68 - 2.04eV for 
t ~ 3eV as shown in Fig. 1. jf represents the FM ex- are given as 
change interaction between LUMO and j'-th tt orbital. 
By using the wave functions obtained in the 7r-electron ■'1(1) = 
approximation, the magnitude of Jf 1 = jf 1 = Jl and 
j¥ = j¥ = Jh are evaluated as 



cesses: J^-d = Ji(i) + Jx{2) + h + J3, and J' 7! _ d = 
2(^i(i) +Ji(2))- These processes are shown in Fig. 4(A) ~ 
Fig. 4(D), and the corresponding exchange interactions 



Jl = [ [ dr 1 dr 2 $J(n)$*(r 2 ) 1 — r #i(r 2 )* 3 (ri), 

J J |ri - r 2 \ 



| 2 |cf \*U p *0MU p , 

3 



and, 



j H ^j2\cpr\ c f^u, 



,, _ 0.01C/ P , 



(11) 



where the coefficient, c*, is the amplitude of the wave 
function of j-th N and C atoms of i-th ir orbital in eq.(l), 
and we have approximated as Un = Uc = U p where Un 
and Uc are intra-Coulomb interactions of 7r orbital on 
N sites and on C sites. We have neglected other interac- 
tions. Then, we find that Jl — AJh- 

Figure 3 shows a schematic picture of the Hamiltonian, 
eq. (2). Since it is assumed that Fe is trivalent (Fe 3+ ), 
and Pc is monovalent (Pc 1_ ) in this Hamiltonian, total 
number of electrons is 5(d orbitals) + l(LUMO) + 4 
(HOMO+l,HOMO+2) = 10. 
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Fig. 3. (Color online) Schematic picture of effective model. The 
green dashed lines indicate the transfer integrals between d or- 
bitals and n orbitals. We neglect HOMO, HOMO+3, HOMO+4, 
HOMO+5, HOMO+6, HOMO+7, d 3z2 _ r 2, and d x2 _ y 2. 



We derive the 7r-d interaction in the Hamiltonian eq. 
(2) by the second order perturbation as follows. 

H n _ d = 2J 7r _ d S i ■ (S»» + S") + 2J^_ d S L ■ S=^2) 

where S L , S yz , S zx , and S xy are the spin operators of 
LUMO, d yz , d zx , and d xy orbitals, respectively. The ex- 
change interaction J 7! - d and J^_ d contain various pro- 
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(A L+1 -A d -U d - 3U' d + j d y - J L 
Jl (tu + *ia) 
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J2 = - 
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(A L+1 -A d - AU' d + 2 J d ) 2 - J 2 L ■ 

2Jh(*?5 + tie) 

72 



(15) 



(16) 



(A d + U d + AU' d - 2J d ) .,„ 

Note that Jin) is AFM, while Ji(2), J2, and J3 are FM 
In the following, we explain the second order perturba- 
tion processes in detail for Ji(i), Ji(2), J2, and J3, re- 
spectively. 



(A) AFM coupling "double Hund's exchange" (J, (1) process) 
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Fig. 4. (Color online) Second order perturbation processes of (A) 
an antiferromagnetic (AFM) interaction process) and (B) 

~ (D) a ferromagnetic (FM) interaction (J^) ~ ^3 processes). 
The dotted lines indicate the hybridization between d orbital and 
7r orbital. 



In Fig. 4, the initial state is (xy) 2 (yx) 2 (zx) 1 for d elec- 
trons. At first, the schematic picture in Fig. 4(A) (the 
J i(\) process ) corresponds to the following process: 
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(1) Initial state: (xy) 2 (yz) 2 (zx) 1 . 

(2) Intermediate state: The electron with down spin in 
the d yz orbital moves to LUMO+1. Because of the 
Hund's rule coupling between d yz and d zx , this in- 
termediate state is favorable. On the other hand, 
if LUMO has a down spin electron, this interme- 
diate state is also favorable because there is a FM 
exchange coupling between LUMO and LUMO+1. 
In this process, we have an energy denominator 
(A L+ i -A d -U d - 3U' d + J d ) 2 - J\ as shown in 
eq. (14). 

(3) Final state: The electron on LUMO+1 moves back 
to the d yz orbital. As a result, we find that the in- 
teraction between LUMO and d zx is AFM. 

Secondly the schematic picture in Fig. 4(B) (the J\{2) 
process ) is 

(1) Initial state: (xy) 2 (yz) 2 (zx) : . 

(2) Intermediate state: The electron with up spin in the 
d yz orbital moves to LUMO+1. Although there is 
no energy gain of the Hund's rule coupling between 
dy Z and d zx , there exists the FM exchange coupling 
between LUMO and LUMO+1. Then, we have an 
energy denominator (A i+1 — A d — U d — 3U d +2J d ) 2 — 
J\ as shown in cq. (15). 

(3) Final state: The electron on LUMO+1 moves back 
to the dyz orbital. As a result, we find that the in- 
teraction between LUMO and d zx is FM. 

Thirdly, the schematic picture in Fig. 4(C) (the J 2 
process ) is 

(1) Initial state: (xy) 2 (yz) 2 (zx) x . 

(2) Intermediate state: The electron with up spin in the 
d zx orbital moves to LUMO+1. Although there is no 
energy gain of the Hund's rule coupling between d yz 
and d zx as in the case of (B), there exists the FM 
exchange coupling between LUMO and LUMO+1. 
In this case, the d yz orbital is the double occupancy. 
Then we have an energy denominator (A^ + i — A^ — 
W' d + 2J d f J 2 . 

(3) Final state: The electron on LUMO+1 moves back 
to the d zx orbital. As a result, we find that the in- 
teraction between LUMO and d zx is FM. 

Finally, the schematic picture shown in Fig. 4(D) (the J3 
process ) is 

(1) Initial state: (xy) 2 (yz) 2 (zx) x . 

(2) Intermediate state: The electron with down spin in 
HOMO+1 moves to the d zx . Because of the FM ex- 
change coupling between tt orbitals, the energy de- 
nominator is (A d + U d + AU' d - 2J d ) 2 - J%. 

(3) Final state: The electron in the d zx orbital moves 
back to HOMO+1. Consequently, we find that the 
interaction between LUMO and d zx is FM. 

From these results, it is found that the process of Fig. 
4(A) is AFM, while the processes of Fig. 4(B) ~ Fig. 
4(D) are FM. 

These mechanisms for exchange interaction are basi- 
cally the Kanamori-Goodcnough (KG) mechanism, be- 
cause they are between two orthogonal orbitals and they 
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use ferromagnetic exchange interactions between 7r or- 
bitals in the intermediate state. 18 However, they are 
much more complicated than the conventional KG mech- 
anism discussed in the TM-oxides, and they can be called 
as "generalized" KG mechanisms. In particular, the pro- 
cess (A) uses two kinds of Hund's rule couplings both 
on the 7r-orbital and on the Fe d-orbital. This process 
leads to the antiferromagnetic coupling Ji(i), in contrast 
to the ferromagnetic one in the usual KG mechanism. 
Thus, we call this process (A) as a " double- Hund's ex- 
change". Since this process is favorable for both of the 
Hund's rule coupling, the energy denominator is gener- 
ally small among the processes in Fig. 4(A) ~ Fig. 4(C). 
As a result, this "double Hund's exchange" (the process 
of Fig. 4(A)) can be dominant in in Fig. 4(A) ~ Fig. 4(C). 
Although it has been believed that the 7r-d interaction 
is caused by a simple super-exchange type interaction, 
we found that the mechanism of 7r-d interaction in Pc 
is a new kind of mechanism because of the absence of 
LUMO-d hopping. 

The effective Hamiltonian in the second order pertur- 
bations for (xy) 2 (yx) 2 (zx) x and (xy) 2 (yx) 1 (zx) 2 , are the 
same, because the amplitude of the transfer integrals be- 
tween d zx and tt orbitals are the same as those between 
dy Z and 7r orbitals. On the other hand, when the initial 
state is (xy) 1 (yz) 2 (zx) 2 , there exists only the process as 
in (A) and (B) of Fig. 4. Then, the exchange interaction 
is given by J'^_ d . 

Figure 5(a) shows the absolute values of Ji(i), j 
J2, and J3 (eqs.(14) ~ (16)) as a function of A d /t d N for a 
parameter set U d /tdN — 5.0, U d /t&$ = 3.0, J d /t<m = 1-0 
and Ji/tdN = 0.1. It is noted that magnitude of var- 
ious matrix elements associated with Coulomb interac- 
tions may be estimated quantitatively by the quantum 
chemical calculations. However, since these parameters 
of Fe-Pc system has not been studied in detail, we as- 
sume a range of parameter values, which we think rea- 
sonable based on the present understanding of transition- 
metal oxides. We can see that J t m is the largest among 
these exchange interactions for A d /t d N > —13, while J3 
is the largest for A^/idN < —13. It is noted that there 
is a restriction of the value of A d /t d N as shown below: 
we estimate the total energies of (3d) 4 , (3d) 5 , and (3d) 6 
state in the zeroth-order Hamiltonian of H pc + H d . When 
A d /t d N < —13.7, the electronic state (3d) 6 becomes 

more stable than (3d) 5 , which is inconsistent with the 
experiments. Similarly, when A d /t d N > —11.7, (3d) 4 be- 
comes more stable than (3d) 5 , inconsistent with the ex- 
periments. Therefore, the parameter region of the (3d) 5 
state is —13.7 < A d /t d M < —11.7 which is shown in 
Fig. 5(a). 

In order to check the validity of the second-order per- 
turbation, we evaluate the energy difference between 
S = and S = 1 by a numerical diagonalization of 
the total Hamiltonian H = H t + H d + Hp c . Figure 5(b) 
shows the result (black dashed-dotted line) as a function 
of A d /t d N- One finds that the ground state is S = 1 
state for —13.7 < Ad/tdN < —12.5, and the interac- 
tion is FM, while the ground state is S — state for 
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Fig. 5. (Color online) (a) A-d/t(a)dN dependence on Ji(i), Ji(2)> 
J2, and J3 at U' d /t = 3.0, Jd/t = 1-0. The blue (red) lines indi- 
cate a ferromagnetic (an antiferromagnetic) interaction. The re- 
alistic parameter region of A^/f^jv is —13.7 < A^/t^jv < —11.7. 
(b)Energy difference between 5 = and 5=1 states as a func- 
tion of Ad/tdN- The black dash-dotted line and the red solid 
line indicate the results of numerical diagonalization and second 
order perturbation, respectively. 



-12.5 < A d /t dN < -11.7, and the interaction is AFM. 
This indicates that the interaction between LUMO and 
d yz (d zx ) orbital can be FM or AFM depending on the 
parameters. For comparison, the energy difference be- 
tween 5 = and S = 1 is also calculated in the second- 
order perturbation by J^-d- This is shown in Fig. 5(b) 
by a red line. It is found that the result of perturbation is 
consistent with the result of the numerical diagonaliza- 
tion. The reason why these results are in good agreement 
with each other is that many body effects in d electrons 
is much stronger than the transfer integrals. We will dis- 
cuss and determine the sign of the exchange interaction 
by comparing with experimental results in §6. 

4. Electronic state of d electrons: Analysis of the 
g-values 

In the previous section, we have discussed the mecha- 
nism of 7r-d interaction. In order to make a detailed com- 
parison with experiments, it is necessary to take account 
of the anisotropy of 7r-d interaction observed experimen- 
tally. This anisotropy is due to the spin-orbit interaction 
and the CEF in the real system. For the discussion of 
the anisotropy, we consider a material, PNP[FePc(CN) 2 ], 
which has a local moment S — 1/2 of Fe ion, while it does 
not have tt electrons (i.e, band insulator). In this mate- 
rial, we can understand the electronic state purely of d 
electrons on Fe site. Thus, in this section we discuss the 



electronic state of d electrons for this PNP-compound 
including its anisotropy, by analyzing the anisotropy of 
the g-values obtained in ESR experiment. 2 

Figure 6(a) shows the schematic picture of the CEF 
splitting of the d orbitals of PNP[FePc(CN) 2 ], where the 
energy splittings are estimated on the basis of the usual 
CEF theory. 19 Since it has been known that Fe 3+ is in 



Pc 2 " Fe 3+ (3d 5 S=l/2) 
Energy _3z^ — 

~0.2eV 



X 2.y2 



LUMO 



= ~0.9eV 

J — . 

zx, yz ;;;;;; 
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-O.leV 



(a) CEF 



(b) Spin-Orbit (c) Magnetic field 
coupling 



Fig. 6. Schematic picture of the energy levels in t2 9 orbitals with 
(a)the crystalline electric field (CEF) effect, (b)the spin-orbit 
coupling , and (c) the magnetic field, comparing with Pc 2— elec- 
tric structure. 



a low spin state from the static magnetic susceptibility 
measurements, five electrons occupy t^g (d xy , d yz , and 
d Z x) orbitals as shown in Fig. 6(a). 

Since the energy difference between i 2 g and e g states 
is about leV, we consider only i 2g orbitals to discuss the 
g-values. Then, an effective Hamiltonian becomes 



-Htot. — Hn 



H p 

magi 



(17) 



where H , H so , and -ff£, ag represent one-body energy 
level, spin-orbit coupling, and the magnetic field along 
the p direction (p = x, y, z), respectively. We write H as 



yz IL yza 



A 



zx ,L zxa\ 1 



(18) 



where A yz and A zx are one-body levels of d yz and d zx 
measured from d xy , respectively. Note that the amplitude 
of A yz is generally different from that of A zx because 
of the distortion in the local structure around Fe (See 
Fig. 7(b)). The spin-orbit coupling in the manifold of ti g 
orbitals is given as 20 



. A 



(19) 



£,m.n 



where a n is the n-th component of the Pauli matrix, 
i.e. <7i = a xi 02 = cry, and 173 — a z . A is a spin-orbit 
coupling constant, and Q mn is the Levi-Civita symbol. 
dt a is the annihilation operator with spin a on the £ 
orbital. Finally, The Hamiltonian of the magnetic field is 



mag mag 



mag' 



(20) 
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where p = 1 ~ 3 (1 = x, 2 = y, and 3 = z). ifmag and 
Hiuag are 

#m*a g = ^ d l d l<r'(<r P U'> ( 21 ) 

^mag = i^P e tmp d\ a d ma , (22) 

for the ^2g orbitals, where go = 2, /i p = [IbH p , and /its 
and -Hp is the Bohr magneton, and the external magnetic 
field along the p axis, respectively. 

When five electrons occupy the t 2g orbitals, there 
are six configurations. On the basis of these configura- 
tions, the eigenvalues and eigenstates of the Hamilto- 
nian eq.(17) are obtained by a simple diagonalization. 
As shown in Fig. 6(b), these eigenvalues consist of three 
two-fold degenerate states (Kramers doublets) without 
magnetic field. Under the magnetic field, the Kramers 
doublet splits as shown in Fig. 6(c), where E+ and E_ 
are the highest energy state and the second highest en- 
ergy state, respectively. Then, the g-values along the p 
(p = 1 ~ 3) axis are given by g p = d[E + — EJ\/ dh p \ h _^ Q . 

Figure 7(a) shows the numerical results of g p (p = 
1 ~ 3) as a function of A yz with A zx being fixed at 
A zx = 0.05 eV and A = O.OleV. At A yz = 0.045eV 



a) b) 




Fig. 7. (Color online) (a) A yz dependence on g Xt g y and g z at 
A zx = 0.05eV and A = O.OleV. The dotted line indicates the 
location of A yz = 0.045eV. (b) Schematic picture of the stacking 
of FePc(CN) 2 unit. 



(the dotted line in Fig. 7(a)), the g-values are evalu- 
ated as g x ~ 0.5, g y ~ 1.3, and g z ~ 3.7, respectively. 
These values are in good agreement with the g-values of 
PNP[FePc(CN) 2 ], g x ~ 0.5,g y ~ 1.1, and g z ~ 3.6, ob- 
tained experimentally. 2 As shown in Fig. 7(b), the differ- 
ence between A yz and A zx will be due to the distortion 
of the apical cyanides caused by the Coulomb interaction 
between the cyanide and the nearest Pc. 

5. Effective model of the 7r-d interaction 

Although the mechanism of ir-d interaction is clarified 
in §3, the sign of the n-d interaction has not been deter- 
mined. In this section, we discuss its sign by calculating 
the susceptibility of a single Fe-Pc system which has one 
localized spin on Fe and on the LUMO, respectively. We 
use an effective Hamiltonian 

H-k-a = Ho + H so + -ff^-d, (23) 



where H , H so , and H^-d are eq. (18), eq. (19), and 
eq. (12), respectively. Note that, in deriving H^-d, we 
have neglected the effects of CEF and spin-orbit coupling 
assuming that these effects gives subsidiary modifications 
of exchange interaction. In eq. (23), we discuss the n-d 
exchange interaction in the presence of Hq and H so , in 
order to make a comparison with experiment. 

As shown in Fig. 8(a), the ground state of d elec- 
trons with the spin-orbit coupling and CEF effect con- 
sist of three Kramers doublets (Same as Fig. 6(b)). We 
express the electron in the highest Kramers doublet as 
a pseudo spin. 21 We introduce the pseudo spin opera- 
tors as S* = jS d , S d = ~/S d , and S* = j'S d , where 

S d = S yz + S zx + S x Vj 7 {X/A zx )l, and i - 1. It 

is noted that the base of pseudo spin S d is the highest 
Kramers doublet. We also find Sf> ~ 7^, S% y ~ 7^, 
and Sf = i'Si where 7" - (X/A zx ) 2 . The absolute val- 
ues of 7 and 7" are of the order of 0.01 for A zx = 0.05eV 
and A = O.OleV. Note that LUMO is the part of the 
conduction band constructed from n orbitals in actual 
materials. Therefore, we also introduce a pseudo-spin op- 
erators of LUMO as S L = aS L where a is a parameter 
reflecting the electronic state of the conduction band. 

By using these pseudo spins, the effective Hamiltonian 
is given by 

H e JJ d ~ 2J 17 (^ + ~SyS d y ) + 2J^ d S^Si (24) 

where J1 = aJ\ = a(Ji(i) + Ji(2)), -h-d = aJn-d, and 
we have set as 7' + 7" = 1. The red wiggly line of Fig. 
8(a) indicates the interaction between S L and S d (n-d 
interaction). 

The effective Hamiltonian of the magnetic field is given 

by 

H z mag = \gii*B<*St + g d zf i B S d z ]h z , (25) 

where h z is the magnetic field along z axis, and g z and 
g z are the g-values of pseudo spin on the ir orbitals and 
d orbitals, respectively. 

Figure 8(b) shows the temperature T-dependence of 
the susceptibility by a simple diagonalization for a pa- 
rameter set Jn-d = 100K (the AFM interaction), J\ = 
200i^, 7 = -0.05 and g% = 2.0 with a = 0.2, 0.3, and 
0.5. Here, we set g factor as g x ~ 0.5, g y ~ 1.1, and g z ~ 
3.6 from the experimental result of PNP[FePc(CN)2]2- 
For T > 5K, the T-dependence is Curie's law, while 
for T < 5K the susceptibility is constant irrespective 
of a due to the antiferromagnetic coupling. The inset 
of Fig 8(b) shows the T-dependence for a parameter 
set J^-d = -2Q0K (the FM interaction), J 1 = lOOiT, 
7 = -0.05, g d z = 3.6 and g% = 2.0 with a = 0.5. The 
amplitude of FM interaction is taken as the same order 
of the result of Ref. 7 The T-dependence is Curie's law 
over the whole temperature range, because the ground 
state is a doublet of = (1, ±1) where S* and S* 

are a total pseudo spin and z component of the pseudo 
spin. 

Here, let us compare the above theoretical result 
with experiments. TPP [Feo.o7Coo.93Pc(CN) 2 ]2 is a suit- 
able system for comparison, because it contains di- 
lute Fe local spins. It is reported that the suscepti- 
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Fig. 8. (Color online) (a) Schematic picture of the n-d interaction 
between the pseudo spin on LUMO (S L ) and the pseudo spin 
on d orbitals (S d ). a is the parameter reflecting the electronic 
state of the conduction band, (b) Temperature dependence of the 
susceptibility in the case of the antiferromagnctic 7r-d interaction. 
Inset: Temperature dependence of the susceptibility in the case 
of the ferromagnetic 7r-d interaction. 



bility of TPP [Fe .o7Coo.93Pc(CN) 2 ]2 increases as T 
decreases, and then it decreases in the range T < 
10K. 5 Since the T-dependence of the susceptibility of 
TPP[CoPc(CN) 2 ]2 is small and almost constant except 
for the low temperature, 4 the origin of the T-dependence 
of TPP[Fe .o7Coo.93Pc(CN) 2 ]2 must be due to the Fe lo- 
cal spin with the 7r-d interaction. By comparing this ex- 
perimental result with the theoretical results shown in 
Fig. 8(b), it is found that the exchange interaction should 
be AFM, because the susceptibility does not diverge at 
the low temperatures. For a ~ 0.3, the T-dependence in 
Fig. 8(b) is similar to the experimental result. 

When we compare the experimental results more 
closely, we find that \ decreases slightly below T ~ 10 K. 
This behavior is not obtained theoretically in Fig. 8(b). 
We expect that this decrease of x is due to the Kondo ef- 
fect which has not been taken into account in our model. 
In order to take account of the Kondo effect, the Fermi 
surface of the 7r-band should be included, which remains 
as a future problem. In this case, the charge dispro- 
portionation (4k F CDW) observed in TPP[CoPc(CN) 2 ] 2 
should be also taken into account. 4 Although these con- 
siderations are necessary for the detailed comparison, the 
conclusion of the AFM exchange interaction between d 
electrons and it electrons will not be changed. 

As discussed above, we found that the 7r-d interaction 
is AFM from the experimental result, while it has been 
suggested that the 7r-d interaction is FM in the quantum 
chemical calculation. 7 Since the AFM interaction occurs 
by the second-order perturbation using the high energy 
states (LUMO+1 and LUMO+2), we suspect that these 
excited states are not taken into account enough in the 
quantum chemical calculation. The extended Kondo lat- 
tice model assuming FM or AFM 7r-d interaction has 
been suggested to be the origin of the negative mag- 
netoresistance. 6,8,9 From our study, one finds that the 
AFM Kondo lattice model suggested by ref. 6 is in agree- 
ment with our result. However, because an isotropic ir- 
d interaction is assumed in ref., 6 it is needed that the 



anisotropy of d electrons in TM is considered. 

6. Conclusion 

In this paper, we have clarified that the mechanism 
of 7r-d interaction in Pc is an extended version of the 
Kanamori-Goodenough mechanism including the double- 
Hund's exchange mechanism that is antiferromagnetic. 
This is different from the super-exchange mechanism oc- 
curring from the hybridization between two orbitals. We 
found that the exchange interaction between d orbitals 
of Fe and LUMO in Pc can be either the FM or AFM 
interaction depending on the parameters on the basis of 
the second order perturbation and numerical diagonal- 
ization calculations. By the comparison with the theo- 
retical result of susceptibility on the basis of the effective 
model with the anisotropy of exchange interactions, we 
concluded that the exchange interaction is AFM. 
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